Reconstructing Dark Energy : A Comparison of Cosmological Parameters 
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A large number of cosmological parameters have been suggested for obtaining information on 
the nature of dark energy. In this work, we study the efficacy of these different parameters in 
discriminating theoretical models of dark energy, using both currently available supernova (SNe) 
data, and simulations of future observations. We find that the current data does not put strong 
constraints on the nature of dark energy, irrespective of the cosmological parameter used. For future 
data, we find that the although deceleration parameter can accurately reconstruct some dark energy 
models, it is unable to discriminate between different models of dark energy, therefore limiting its 
usefulness. Physical parameters such as the equation of state of dark energy, or the dark energy 
density do a good job of both reconstruction and discrimination if the matter density is known to 
high accuracy. However, uncertainty in matter density reduces the efScacy of these parameters. A 
recently proposed parameter, Om, constructed from the first derivative of the SNe data, works very 
well in discriminating different theoretical models of dark energy, and has the added advantage of 
not being dependent on the value of matter density. Thus we find that a cosmological parameter 
constructed from the first derivative of the data, for which the theoretical models of dark energy 
are sufficiently distant from each other, and which is independent of the matter density, performs 
the best in reconstructing dark energy from SNe data. 

PACS numbers; 



I. INTRODUCTION 

The nature of dark energy is one of the most tantalizing 
mysteries in cosmology today. Observations of high red- 
shift type la supernovae tell us that the expansion of the 
universe is accelerating at present which can not be 
satisfactorily explained in the standard cold dark matter 
(CDM) scenario. Independent observations of the cos- 
mic microwave background 2] and large scale structure 
tell us that two-thirds of the present density of the 
universe is composed of some unknown component. The 
study of this unknown "dark energy" is of great interest 
among cosmologists today. 

Various theoretical models for dark energy have been 
suggested, the simplest being the cosmological con- 
stant model with constant dark energy density and equa- 
tion of state wde — Other models of dark energy 
include physically motivated models like the scalar field 
quintessence and Chaplygin gas models, as well as ge- 
ometrically motivated models like scalar-tensor theories 
and higher dimensional braneworld models (see [3l and 
references therein). 

The growing number of theoretical dark energy mod- 
els has inspired a complementary, data-driven approach 
in which the properties of dark energy are reconstructed 
from the data by studying the cosmological parameters 
characterizing dark energy. Two primary methods are 
used for reconstructing cosmological parameters. In the 
first approach, known as parametric reconstruction, a 
sufficiently general fitting function is used to represent 
the parameter in the analysis. This suffers from the pos- 
sibility of bias, depending on the form chosen for the pa- 
rameter. The second method is that of non-parametric 
reconstruction, in which no specific form is assumed for 



the parameter. The difficulty with this is that the param- 
eters of interest are usually obtained by taking the first 
or second derivative of the data, therefore, direct recon- 
struction involving differentiation of noisy data can lead 
to large errors. Many different cosmological parameters 
have been suggested both these reconstruction methods 
(see JS] and references therein). In this work, we attempt 
to study the relative efficacy of the different cosmologi- 
cal parameters in reconstructing dark energy from obser- 
vations, and discriminating between different theoretical 
dark energy models, using the parametric reconstruction 
formalism. 

The paper is arranged as follows- section |lT] contains a 
description of the data and methods used in the analysis, 
section IIIII outlines the results, and section IIVI presents 
the conclusions. 



II. METHODOLOGY 



This work attempts to classify the different cosmologi- 
cal parameters that characterize dark energy in terms of 
their efficiency in constraining the nature of dark energy. 
In order to do this, we analyze cosmological observations 
using the different parameters to compare how accurately 
these parameters reconstruct and discriminate between 
various theoretical models of dark energy. In this analy- 
sis, we primarily use Type la supernova data along with 
information on the present day matter content of the uni- 
verse. 



2 



A. Supernova Data 

Type la supernova are the most direct evidence for the 
existence of dark energy at present. From early twenti- 
eth century, they were investigated as standard candles 
and many attempts were made to use them to measure 
the Hubble parameter and the deceleration of the uni- 
verse [6]. The first cosmologically significant results for 
deceleration of the universe were produced in the late 
nineties, when two observational groups 0, @ indepen- 
dently showed that the expansion of the universe was ac- 
celerating. Since then, there have been numerous other 
SNe surveys P, Q, and despite being plagued by sys- 
tematics, these remain our best observational tool for 
studying dark energy. 

Supernova data by itself is insufficient to break the de- 
generacy between dark energy parameters and the curva- 
ture of the universe. Since the objective of this paper is 
to constrain dark energy using SNe data, we restrict our 
analysis to a flat model (17^ = 1) of the universe which 
is favoured by the current CMB data [2]. The data is in 
the form 



^b(z) = 51ogiodL(z) + M 



(1) 



where A4 represents a noise parameter usually marginal- 
ized over, and the luminosity distance dL{z) is related to 
the cosmological parameters in a flat universe as- 



dL{z) = c(l + z) 



dz 
H{z) 



H{z) = - = i/oVfWT + ^F+^W^ 

a 



(2) 
(3) 



(a is the scale factor representing the expansion of the 
universe, H{z) is the Hubble parameter, c denotes the 
speed of light and Hq the present value of the Hubble 
parameter in km/s/Mpc). 

We perform a maximum likelihood analysis on the su- 
pernova data to obtain constraints on the various dark 
energy parameters, the likelihood being defined as 



C oc e-^''"^ 



X 



E 



(4) 
, (5) 



where pj are the parameters of the fitting function chosen 
to represent the cosmological parameter being studied. 

We use one of the most current SNe datasets for our 
analysis. However, as we shall see, the current data is not 
yet of such a quality that it could strongly discriminate 
between different models of dark energy. We therefore 
also simulate three datasets based on our expectations 
from future surveys p!o| . to study the information that 
could be obtained from future data. The datasets used 
in our analysis are- 

• Dataset I : Currently available Union2 dataset, 
with ~ 550 SNe between redshifts z — — lA, and 
average statistical error of cr^^ ^ 0.1 — 0.3 mags. 



• Dataset II A : Simulated dataset based on future 
JDEM-like SNe surveys containing - 2000 SNe dis- 
tributed over a redshift range of z = — 1.7 with a 
larger concentration of supernovae in the midrange 
redshift bins (z ~ 0.4 — 1.1) and average statisti- 
cal errors of a^^^ = 0.13 mags [l3|. The theoretical 
model used to simulate the data is the cosmological 
constant model with Hq = 72 km/s/Mpc, the mat- 
ter density fiom = 0.27, and the equation of state 
of dark energy wde = ^ 1- 

• Dataset II B : Simulated dataset based on a JDEM- 
like survey as in II A, using a theoretical model of 
quintessence with a minimally coupled scalar field 
whose equation of motion is given by 



dV 



= 



(6) 



with a potential [ll| 

F(</)) = Fo(/." ; a = 2 , (7) 
and the equation of state of dark energy 



Wde 



(8) 



Ooni and -ffo a-re the same as in Dataset II A. 



• Dataset II C : Simulated JDEM-like dataset using 
a variable dark energy model with the equation of 
state given by [l2| 



Wde{z) = Wo + {Wm-Wo)) 



1 + e^t(i+-t) 
1-e^ 



g37 _|_ g At(l + 2t) 



e At(i+^) ^ gAt(i+it) 



(9) 



with the values wq = —1.0, Wm — —0.5, zt = 
0.5, At = 0.05. This model has wde > -1 
everywhere, and may be realized by a standard 
quintessence field, ilom and Hq are the same as 
in Dataset II A. 



B. Parameters of interest 

Several parameters have been suggested in the litera- 
ture for reconstructing the nature of dark energy. These 
can be broadly divided into two categories- 

• Geometrical parameters of dark energy: These are 
parameters that can be constructed directly from 
the scale factor a and its time derivatives. Exam- 
ples are the Hubble parameter H{z), and the quan- 
tity Om(z) [l3|, both constructed from the first 



3 



derivative of the supernova data- 
H{z) 



a 
a 



0m(2) = 



(1 + zf - 1 



(10) 
(11) 



the deceleration parameter q{z), constructed from 
the second derivative of the data- 



q{z) 



-1 



dlogiJ 
(ilog(l + z) 



(12) 



and the Statefinder parameters [14| . constructed 
from the third derivative of the data- 



r{z) = 
s{z) = 



a a 



a"' 

r - 1 
3(^ 



(13) 
(14) 



• Physical parameters of dark energy: These pa- 
rameters, in addition to the scale factor and its 
derivatives, also contain physical information such 
as the matter density. Examples are the dark en- 
ergy density normalized to the critical energy den- 
sity poc — SHq/SttG, which we denote as iloEiz), 
constructed from the first derivative of the data- 



^^DEiz) 



IP 

Hi 



nOn,{l+zf 



(15) 



and the equation of state of dark energy wde(z), 
constructed from the second derivative of the data- 



wde{z) 



DE 



2(l + z)/3 (din H/dz)~l 
1 ~ {Ho /Hy non.il + z)3 



(16) 



The disadvantage for these parameters is their 
dependence on the physical nature of dark en- 
ergy, which means that for non-physical models of 
dark energy which do not follow the standard Ein- 
stein equations (such as modified gravity models, 
or DM-I-DE interacting models), these parameters 
may not be very well-defined. Also, these param- 
eters depend on the matter density, and therefore 
their reconstruction may have added bias due to 
uncertainty in our knowledge of the matter density. 

To study the efficacy of the different parameters in 
quantifying dark energy, we use both geometrical and 
physical parameters- we select the geometrical param- 
eters Om(z) and q(z), and the physical parameters 
^deIz) and wde{z)- For each parameter we need to 
choose a fitting function which will be used in the like- 
lihood analysis. We experimented with different fitting 
functions, e.g., comparing the polynomial expansion in 
redshift z with that in scale factor a for wde{z); and in 



each case we choose the fitting function which results in 
the best reconstruction of the parameter of interest. The 
fitting functions chosen for the different parameters are- 



Oni(z) = Omo + 0mi(l + z) 

. ^ , qiz 

cliz) - 'Zo + 

^De{z) = (1 - f^Om - f^l - f^2) 

+ Z) + VL2{1 + zf 

= 0.0 + Vliz + Vt2Z^ 

WlZ 



wde{z) = Wo + 



l + z 



(17) 
(18) 



(19) 
(20) 



of which the fitting function for q{z) was first introduced 



in pL5|, that for Q,r)E{z) was introduced in [14|, and that 
for wde{z) in [l6|. We use these fitting function for 
the likelihood estimation to obtain confidence levels on 
the cosmological parameters, which can then be used to 
discriminate between different theoretical models of dark 
energy. 

For the geometrical parameters of dark energy, no fur- 
ther information is necessary for performing the analy- 
sis. However, for the physical parameters fluE{z) and 
wde{z), information is required on the matter density. 
We therefore do the analysis for the physical parameters 
by first fixing the value of matter density to r^om = 0.27 
(which is the true value of ilom for the datasets II A, 
B, C, and is the expected value today from large scale 
structure), then by marginalizing over the matter den- 
sity using riom = 0.27 ± 0.03 Q. The results for the 
second case are expected to be worse for the physical pa- 
rameters, while for the geometrical parameters there is 
no difference since they do not depend on fiom- 



III. RESULTS 

We first study the currently available data. The results 
for the Union2 dataset, i.e., dataset I, using riom — 0.27, 
are shown in figurc[TJ We see that the 2a confidence levels 
for all four parameters are extremely large, and the three 
very different dark energy models proposed in datasets 
II A, B, C all fall within these confidence levels for at 
least a large part of the redshift range. Therefore, with 
the current data, at 95% confidence, three widely varying 
models of dark energy cannot be distinguished very well 
from each other using different cosmological parameters, 
even if the matter density were known exactly. 

Since the current data cannot be used to effectively 
constrain theoretical models of dark energy, we now look 
at the constraints that may be obtained from future 
observations. Dataset II A is a JDEM-like simulated 
dataset for which the true underlying model is the cos- 
mological constant. Figure [5] shows the 2a confidence 
levels on the different parameters for this dataset. We 
see that the parameter Om has the lowest errors and 
represents the true model quite accurately. The deceler- 
ation parameter q also has narrow errors, however, since 
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FIG. 1: Variation of different cosmological parameters with redshift for dataset I (real data). Panel (a) shows the quantity 
Om{z), panel (b) shows ^1de{z), panel (c) shows q{z) and panel (d) shows wde{z). The thick solid lines represent the best-fit, 
and the solid grey contours represent the 2a confidence levels. The dashed, dotted and dot-dashed lines represent the true 
model for datasets II A, B and C respectively. In panels (b) and (d), the solid grey contours represent the 2a confidence levels 
when Jlom is known exactly. 



the three cosmological models considered are much closer 
in g-space than in the other parameter spaces, even with 
these narrow errors, q cannot distinguish between the 
dark energy models of dataset II A and dataset II B. 
Both the equation of state wde and the energy density 
floE have somewhat larger errors but do a better job as 
compared to the deceleration parameter of constraining 
the true model as well as distinguishing it from the other 
two dark energy models considered, when exact value of 
ilom is known. The parameter Qde has especially nar- 
row error bars at low redshifts, but after a redshift of 
z ^ 0.7 the error bars become larger. When value of 
Oom is not known exactly this leads to additional uncer- 
tainty in the results, as well as the possibility of bias. A 
higher value of flom achieves a luminosity distance sim- 
ilar to that obtained by dark energy with an increasing 
w, therefore there is a degeneracy between riom and the 
dark energy parameters, leading to the larger error bars 
on the parameters wde and floE- 

Next, we look at dataset II B, which has an underlying 



model of quintessence with a slowly varying equation of 
state. Figure [3] shows the 2a confidence levels for the dif- 
ferent parameters for this dataset. The results are very 
similar to that for dataset II A. The parameter Om recon- 
structs quite accurately with narrow error bars at high 
redshifts, but is slightly less effective at z = 0. The pa- 
rameter q reconstructs accurately at low z at less so at 
high z. The parameter Qde is excellent at reconstruction 
at low redshift when ilom is known, while w reconstructs 
reasonably as well when fJom is known. When ilom is not 
known accurately, the reconstruction has a bias for floE 
and w, and the error bars are larger. 

Dataset II C uses a model with a stronger variation 
of the equation of state, from wde = ^1 today to 
Wde = —0.5 at z ~ 1. The results for this dataset 
are shown in figure S) We see that both the deceler- 
ation parameter q and the equation of state wde are 
unable to reconstruct the true model (dot-dashed lines) 
accurately, because the parameterization chosen is not 
general enough to reconstruct the true model accurately. 



(a) (b) 
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FIG. 2: Variation of difFerent cosmological parameters with redshift for dataset II A. Panel (a) shows the quantity Om(z), 
panel (b) shows Qde{z), panel (c) shows q{z) and panel (d) shows wde{z). The shaded contours represent the 2a confidence 
levels, the thick solid lines represent the best-fit, the thick dashed lines in panels (b) and (d) represent the best-fit when Qom is 
marginalized over . The dashed, dotted and dot-dashed lines represent the true model for datasets II A, B and C respectively. 
In panels (b) and (d), the grey solid contours represent the 2a confidence levels when Qom is known exactly, while the hatched 
contours are marginalized over Qom ~ 0.27 ± 0.03. 



Even when figm is known, the equation of state appears 
to favour wde < — 1 at z = at 2a, which would lead 
to an extremely different model of dark energy violating 
the Weak Energy Condition. The dark energy density 
i^DE does better in reconstructing the true model when 
51om is known, especially at low redshift. This is because 
the dark energy density, which is obtained from the first 
derivative of the SNe data, has usually less sharp varia- 
tions than the equation of state or deceleration parame- 
ter, which are related to the second derivative of the data. 
However, when ilom is not known exactly, i^DE performs 
much worse, especially at low redshift, because of the de- 
generacy between floE and the matter density. As in the 
previous cases, Om appears to capture the behaviour of 
this model very well. 

As we see from the figures [2l [3l HI the different pa- 
rameters perform very differently in reconstructing the 
models chosen. We quantify the efficiency of the differ- 
ent parameters using a distance criterion which calculates 



the distance of the true model from the best-fit, weighted 
by the 2a errors obtained by the fitting procedure. The 
lower this value, the more accurate the reconstruction, 
and the narrower the confidence levels on the parame- 
ter. A large value of this parameter could mean either 
very large errors, which lowers the discriminatory power 
of the parameter, or a biased fit that is far away from the 
true result. We also quantify the distance criterion for 
the other two models in each case and take the ratio of 
this quantity for the true model to that for the nearest 
of the other models, since this highlights the discrimi- 
natory power of the parameter. The lower this value, 
the better the parameter at discriminating other models 
from the true model. We note that this quantity is of 
course dependent on the models chosen, however, for the 
three models studied here, it can give a comparative as- 
sessment of the discriminatory power of the parameters 
considered. For a different set of models, the value of 
this quantity would obviously change, but the compari- 




FIG. 3: Variation of different cosmological parameters with redshift for dataset II B. Panel (a) shows the quantity Om(z), 
panel (b) shows ^1de{z), panel (c) shows q(z) and panel (d) shows wde(z). The shaded contours represent the 2cr confidence 
levels, the thick solid lines represent the best-fit, the thick dashed lines in panels (b) and (d) represent the best-fit when ilom is 
marginalized over . The dashed, dotted and dot-dashed lines represent the true model for datasets II A, B and C respectively. 
In panels (b) and (d), the grey solid contours represent the 2o confidence levels when f2om is known exactly, while the hatched 
contours are marginalized over fiom = 0.27 ± 0.03. 



son between the cosmological parameters would be sim- 
ilar, especially since the three models chosen are quite 
different from each other in behaviour. For a parameter 
p, the efficiency criterion is defined as- 

= ^^^Mnc - Pi,f,tf(rlf^,^fit) , (21) 

i 

and the discrimination ratio is defined as- 

^^(True Model) 
(52 (Nearest Model) ' ^ ' 

where fr^^^ j^^-j is the 2a error on the reconstructed param- 
eter. 

The above quantities are calculated in Table U for the 
different parameters for all three models. We see that 
for the different parameters, the on the data is quite 
similar for all cases, showing that they all reconstruct suc- 
cessfully in the data space. However, since the parameter 



space is linked to the data space by one or more differ- 
entiations, the reconstruction of the actual parameters, 
quantified by 5"^ (1^ . is quite different for the different 
parameters. The quantity Om performs consistently well 
in reconstructing the dark energy models considered, and 
has the lowest values of 5^ . The deceleration parameter q 
performs well for Models A, B, but does not do so well for 
Model C which is a model not well characterized by the 
parameterization chosen for it. The dark energy density 
and the dark energy equation of state both work reason- 
ably when the matter density is known, however, their 
performance degenerates when uncertainty is introduced 
in riom, especially in Model C which is the most strongly 
varying. The discrimination criterion (last column of ta- 
ble) shows how well the parameter can discriminate the 
true model from the other two models considered. We 
see that the parameter Om performs very well on this 
criterion as well. The deceleration parameter, which per- 
formed reasonably well on J^, does not perform as well on 
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FIG. 4: Variation of different cosrnological parameters with redshift for dataset II C. Panel (a) shows the quantity Om{z), 
panel (b) shows 0,de{z), panel (c) shows q{z) and panel (d) shows wde{z). The shaded contours represent the 2a confidence 
levels, the thick solid lines represent the best-fit, the thick dashed lines in panels (b) and (d) represent the best-fit when f2om is 
marginalized over . The dashed, dotted and dot-dashed linos represent the true model for datascts II A, B and C respectively. 
In panels (b) and (d), the grey solid contours represent the 2a confidence levels when fiom is known exactly, while the hatched 
contours are marginalized over fiom = 0.27 ± 0.03. 



this criterion, because the different cosmological models 
are quite close together when represented by this quan- 
tity. Thus the deceleration parameter is not very useful 
when it comes to choosing one dark energy model over 
the others. The dark energy density ^de and the equa- 
tion of state wde both perform reasonably when flom is 
known exactly, but once again their performance degen- 
erates when priors are put on fiom, since the error bars 
increase. Overall, in both the criteria considered, Om 
appears to perform the best, while ^de performs well at 
low z provided fiom is known, and w performs moderately 
well when fiom is known. 



IV. CONCLUSIONS 

We reconstruct four different cosmological parameters 

to examine their potential in extracting information from 
observations about dark energy. We find that the pa- 



rameter Om, which is constructed from the Hubble pa- 
rameter, has the narrowest errors, and gives the most 
accurate reconstruction. We also see that results for the 
physical parameters Ude, wjj e are extremely dependent 
on knowledge of the matter density, and without very 
tight bounds on matter density, physical parameters of 
dark energy may not perform well in reconstruction dark 
energy. The deceleration parameter reconstructs quite 
well, but has low discriminatory powers, since different 
models of dark energy have very similar values for this 
parameter. 

Overall it appears that the parameters constructed out 
of the first derivative of the data, i.e. the Hubble pa- 
rameter, are somewhat better poised to give information 
about the nature of dark energy than the second deriva- 
tives, and geometric parameters of dark energy, such as 
Om, are better at reconstructing dark energy since they 
arc not biased by lack of knowledge about the matter 
density, provided they are constructed in such a way that 
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TABLE I: Xmin ^^^^ different parameter reeonstructions using the three simulated datascts, the efHeieney eriterion ( |21| > which defines how 
well the model reconstructs the truth, and the discrimination criterion l |22t . which defines how well the model discriminates from the other two 
truths considered. 



Dataset 




Parameterization 


Amm 




5p{Truc Model) 

<52 (Nearest Model) 
E 


II A 


0.27 


Om{z) = Omo ± Orni(l ± z) 2349.48 
q{z) = qo - qiz/(l + z) 2349.77 
QoEiz) = flo + niz + n2Z^ 2349.45 
w{z) = wo-Wiz/{l + z) 2349.44 


lO"'^ 

0.011 
0.019 
0.014 


0.001 
0.668 
0.009 
0.010 




0.27 ±0.03 


flDE{z) = 0,0 + OiZ + Q.2Z^ 
w{z) = Wq — WlZ / {1 -\- z) 


2349.03 
2349.26 


0.164 
0.173 


0.882 
1.395 


II B 


0.27 


Om(z) = Omo ± Omi(l ± z] 
q{z) = go - qiz/(l + z) 

Q.De{z) = Oo ± OiZ + Q.2Z^ 

w{z) = mo — Wiz/{1 ± z) 


2349.75 
2350.08 
2349.54 
2349.45 


0.0002 
0.047 
0.119 
0.042 


0.007 
0.652 
0.010 
0.009 




0.27 ±0.03 


Ode{z) = S7o + Oiz + Q.2Z^ 
w{z) = Wo — Wiz/{1 + z) 


2349.18 
2349.37 


0.292 
0.201 


1.002 
1.292 


II C 


0.27 


Om(z) = Omo ± Omi(l ± z) 2349.75 
q{z) = qo - qiz / [1 + z) 2349.50 
0.de{z) = Oo + 0-iz^Q.2Z^ 2349.46 
w{z) = wo-wiz/{l + z) 2349.97 


0.0001 
0.104 
0.039 
0.077 


0.006 
0.479 
0.007 
0.023 




0.27 ±0.03 


VIde{z) = 0,0 + OlZ ± 0.2Z^ 

w{z) = Wo — Wiz/{1 + z) 


2349.19 
2348.96 


1.981 
1.163 


0.452 
0.356 



they can discriminate between different models of dark 
energy. To obtain information about dark energy from 
the physical parameters of dark energy, it is important 
to have independent sources of observation of the matter 
density. Thus, for a parameter to obtain maximum infor- 
mation out of the data, it works better if it is constructed 
out of the first derivative of the data or less, does not de- 
pend on other physical parameters such as the matter 
density, and is constructed such that different dark en- 
ergy models can be easily discriminated from each other 
with the parameter. 

We note here that the reconstruction methods con- 
sidered here are that of parametric reconstruction, non- 
parametric reconstruction of cosmological parameters 
would suffer less from biases in reconstruction such as 
those Model C. However, typically these methods have 
higher errors, and also the issues of degeneracy with 



other, non-dark energy parameters exist for these meth- 
ods as well. Therefore if a simple cosmological parame- 
ter, like Om, can be constructed which performs well for 
a large class of models, both in reconstructing and in dis- 
criminating between models, and is also independent of 
other non-dark energy parameters such as rJom, it would 
be extremely useful for understanding the nature of dark 
energy. 
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